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Abstract 

We propose an approach to the problem of low but finite temperature dynamical correlation 
functions in integrable one-dimensional models with a spectral gap. The approach is based on 
the analysis of the leading singularities of the operator matrix elements and is not model specific. 
We discuss only models with well defined asymptotic states. For such models the long time, large 
distance asymptotics of the correlation functions fall into two universality classes. These classes 
differ primarily by whether the behavior of the two-particle S matrix at low momenta is diagonal 
or corresponds to pure reflection. We discuss similarities and differences between our results and 
results obtained by the semi-classical method suggested by Sachdev and Young, Phys. Rev. Lett. 
78, 2220 (1997). 

PACS numbers: 71.10.Pm, 72.80.Sk 
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I. INTRODUCTION 



In this paper we discuss the problem of finite temperature time-dependent correlation 
functions in integrable models. We do not consider here one- dimensional models with gapless 
linear spectra, for which there exists a general solution: the Matsubara n-point correlation 
functions at finite temperature can be obtained from the corresponding T = correlators 
by a conformal transformation. Instead we focus on models with spectral gaps. Finite 
temperature dynamics for that class of models is not yet fully understood. The spectral 
gap, M, provides an energy scale separating the high temperature (T >> M) and low 
temperature (T ^ M) regimes. In the former regime the leading order physics reduces 
to that of a conformal field theory at high temperatures. We thus concentrate on the low 
temperature regime T <C M. In this limit we will argue that the relevant physics is governed 
by the zero momentum limit of the interactions between excitations in the theory. These 
interactions refiect the integrable nature of the models: the N-particle scattering matrix 
is a product of the two-particle ones. It is known that the zero momentum limit of the 
two-particle scattering matrix, S{0), in an integrable system can be either diagonal, as it 
occurs, for instance, in the quantum Ising model, or equal to the permutation operator. The 
latter alternative is realized for the solitons of the sine-Gordon model at general values of 
the coupling constant. This dichotomy suggests the existence of two universality classes for 
the low temperature dynamics. As we will demonstrate in this paper, this expectation is 
met for a class of order parameter-like operators. 

The most advanced methods of evaluating correlation functions for integrable models 
w.h spe*a. gap. a.e based o. the fonn fae.o. teeh,.,ue (see fl fl 3 fo. .views). This 
technique allows one to calculate matrix elements (form factors) of the operators between 
the exact eigenstates entering the Lehmann expansion for the correlation functions. At finite 
temperature this Lehmann expansion for the two-point correlation function of the operator, 
O, has the following form: 

G'^ix^t) = 4Tr(e~^^O(x,t)C(0,0)) 
Z 
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^ En,.„;m,s„ e"^""'" (^^ t)\m, s^) {m, g^|0(O, 0) \n, s^) 

Here the state, \n, s„), denotes a set of n-particles carrying quantum numbers, The final 

expression in the above involves a double sum. The first sum appears as a sum weighted 
by Boltzmann coefficients while the second sum arises from inserting a resolution of the 
identity between the two fields. The two-point correlation function has thus been reduced 
to a question of computing and summing form factors. 

The form factor approach turns out to be highly effective for T = two-point functions. 
This is particularly so in computing spectral functions where at finite energies the sum corre- 
sponding to that above reduces to an expression with a finite number of readily computable 
terms. However, its generalization to multi-point functions, as well as for T 7^ faces dif- 
ficulties caused by singularities that appear in the needed form factors 1]. In this paper 
we propose a method of summation of these singularities that is suited to the extraction of 
the long distance, long time, low temperature asymptotics of the correlation functions. In a 
particular case, our results reproduce the semi-classical formulae obtained in Refs. 
In other instances, we find divergences between our formulae and those derived using the 
semi-classical methodology proposed in Refs. 

The outline of the paper is as follows. In Section II we consider the simplest of the 
integrable models, the Ising model. Using expressions obtained in Ref. 1^ from a lattice 
analysis of this model, we analyze the asymptotics of the spin-spin correlation functions. In 
doing so we are able to cast the spin-spin correlation function in a form nearly identical to 
what would be found using the Lehmann expansion. We take this similarity as instructive: it 
both indicates how to handle singularities that appear in the Lehmann expansion and it tells 
us how to proceed with more general integrable models. Thus in the next section. Section 
III, we generalize our computations to arbitrary integrable models. We find two universality 
classes of behavior. One class of behavior is reminiscent of the Ising model. Here we find no 
discrepancies between our approach and that of the semi-classical approach. For the second 
universality class we find novel scaling behavior that is not reproduced semi-classically. We 
consider the implications of these two universality classes for the physics of the sine-Gordon 
model in Section IV. We conclude our paper with a brief discussion of the results. 
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II. FINITE TEMPERATURE SPIN-SPIN CORRELATION FUNCTIONS IN THE 
ISING MODEL 



To understand how to deal with the singularities that appear in the Lehmann expansion, 
we first turn to the quantum Ising model. In the following section we will extend our 
approach to other integrable models. 

We first review the basics of the model. The Hamiltonian, H, of the Ising model is 
equivalent to a Hamiltonian of noninteracting massive fermions 

n n p 

e{p) = ^J{J -hf + AJhsin\p/2), (2) 

where ~ Y\.j<n '^j disorder operator and F, are fermion annihilation and 

creation operators. The fermions represent walls between domains with different orientations 
of the magnetization and so are to be thought of as solitons. 

Calculations simplify in the continuum limit, when the spectral gap M = \ J — h\ is much 
smaller than the bandwidth ~ J and the spectrum is relativistic e{p) = a/ c^p^ + , = 
Jh. In this case, energy and momentum of a quasi-particle are conveniently parameterized 
by a rapidity, 6, {cp = Msinh^). Then the eigenstates of Hamiltonian Q are labeled by 
sets of rapidities, {Oi}, such that the energy and momentum of the system are equal to 

n n 

E = M'^coshdi, P = c^^M^sinh^i. (3) 

i=l i=l 

Below we set c = 1. Operators and /i^ (we will call them a and /i) have infinitely many 
matrix elements. By Lorentz invariance, the form factors depend on rapidity differences j^: 

(^1, ...en\(y\e[, ...e'J = AM^/^ J]tanh(%/2) J]tanh(e;^/2) J]coth[(e, - ± te)/2], (4) 

i<j p<q i,p 

where aX h < J (the ordered phase) n + m is even for a and odd for yU (for /i > J it is 
the other way around). Here 9ij = 9i — 9j and A is a known numerical constant. The 
singularities in (jl)) appear as soon as some in and out rapidities coincide. Singularities of 
this kind are known as kinematical poles and appear routinely in integrable models. There 
are operators however where the matrix elements are free of singularities. In the Ising model 
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such an operator is the energy density operator, 

a^{x) = Y,e'''"likMk-q)FkF,.k, h = ^{k) = + k/ e{k). (5) 

k 

We will not in general consider the correlation functions of such operators in the paper 
(although see Section IVc). We can say however their long time, large distance behavior 
will generally be ballistic. 



A. Possible Approaches 

In the Ising model the form factors are known explicitly and it has been a common belief 
that their analysis would lead to a better understanding of the general problem. Several 
approaches to the problem of the singularities that appear in the Lehmann expansion are 
possible: 

• One can try to regularize the singularities by considering a finite system. Understand- 
ing the singularities as arising from an infinite sized system has been advocated both 
in Ref. ^3 in Ref. j^. In the Ising case it is clear how this regulation would 
proceed. In the infinite system the singularities arise when the momenta of the in 
and out states coincide. However if the system is treated as on a ring of finite size, 
the in and out states must be quantized differently. The spin operator in the Ising 
model connects the Ramond and Neveu Schwarz sectors of the Hilbert space. In the 
Ramond sector, the fermionic modings are periodic while in the Neveu Schwarz sector 
they are anti-periodic. (These features of the Ising model Hilbert space are reviewed 
in Ref. j^.) Thus at finite system size the momenta simply cannot coincide and 
the expressions are singularity free. However the summation of form factors is then 
discrete not continuous. To the best of our knowledge no one has then succeeded in 
explicitly performing the necessary double summation. 

• In a similar spirited approach, one can interchange imaginary time and space axes 
using a Wick rotation. The finite T Matsubara Green's functions of the infinite system 
then become zero temperature correlators for a system on a circle of circumference, 
1/T. While the corresponding form factor expansion is manifestly free of singularities, 
one faces the problem of calculating form factors and excitation energies in a system 



5 



of finite size. For the Ising model this was achieved in Ref. Q]. The two-point 
correlation functions of cr's and /i's were obtained rigorously in the form of series. 
Below we use these series to extract the long time and distance asymptotics of the 
dynamical correlation functions. 

One can calculate correlations by re-expressing them as the determinant of an integral 
operator of the Fredholm type. In this way results were obtained for the spin-1/2 XX 

n n 

magnet and the Bose gas [13||. For a general exposition on this methodology see 
Ref. Q. 

Exploiting determinantal expressions of the matrix elements arising from explicit ex- 
pressions for the Bethe eigenfunctions, Refs. Q| have obtain expressions for the 
zero temperature correlation functions in the Heisenberg XXZ chain. It is likely this 
approach can be readily extended to finite temperature computations. 

In an approach adopted in Ref. one can imagine treating the computation of cor- 
relation functions at finite temperature in terms of a renormalized 'thermal' vacuum 
and a renormalized set of 'thermal' excitations. The computation of correlation func- 
tions then proceeds in a fashion similar to that of the zero temperature calculation 
but with the form factors obeying a revised set of axioms. More recently, the notion 
of a thermal vacuum was explored in the work of Ref. {h | . 

The sole method that does not trade on the underlying integrability of the models 
of concern was suggested by Sachdev and Young in the form of a "renormalized 
classical" approach valid at finite T and long time. The result for the causal correlation 
function of the spin field in the ordered phase valid for T <^ M is 



D{x,t) = {a{t,x)a{0,0)) =CM^/^exp 



7T dp 



(6) 



where e = a/M^ + ^ M + p'^/2M and C is a numerical constant. This expression 
suggests that the typical correlation length scales as t. One can approximate © as 
D{x,t) ^ Cexp[— 'n,(T)max(|x|, f |t|)] where v = ^Jt^tJ^M is the thermal velocity of 



kinks and n(T) = ^JTMJt^q-^^I^ is the average number of solitons (fermions). For 



the case of the Ising model, we show below that 
form factor approach using the results of Bugrij 



(jni) also follows directly from the 
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• The final possible approach involves attacking directly the double sum Lehmann ex- 
pansion (Eqn. 1) through the use of infinite volume form factors. To understand how 
this is to be done in general, we employ the series expansions from Ref. which will 
provide the necessary clues to understanding the singularity structure in this case. 



B. Long space and time asymptotics of the spin-spin correlation function 

In the section we use the results of Ref. jl^ to understand the asymptotics of the spin- 
spin correlation function in the Ising model. From this we will be able to follow a natural 
path to generalizing these results to integrable models at large. 

The result for the Matsubara Green's function of the Ising model order parameter field 
obtained in Ref. [l^ is 

(a(r,x)a(0,0))= (7) 

Af=l ^ >' qi,...q2N i=l ' 

where r is imaginary time, 

n(T) = j ^e-^^(^) ^ VrM/27re-*^/^, (8) 




is the soliton density, q = 2TiTm (m integer), and e(g) = a/ + q^. The term in the 
exponent, 

Viq) = ^l ^,^^y^lncoth[/3e(x)/2], (9) 

at T ^ M is exponentially small and will be neglected from hereon. The symmetry breaking 
transition at T = leads to a finite mag netization, (a) = ±[CM^/^Y/'^. This is reflected 
in the zeroth order term in Eq. (jZj). We see that this expression for the correlator has the 
form of a zero-temperature Lehmann expansion: the sum over N represents intermediate 
states with N particles while the sum over the g^'s are momentum sums for a given set of 
N particles. The rightmost term in Eq. ((7j) can then be interpreted as a finite volume form 
factor squared. We can see that the correlation function, Eq. ((7j), is periodic in r, as it must 
be. The entire expression can then be thought of as a Wick rotation of a computation done 
in the picture with periodic boundary conditions in the space direction. 
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FIG. 1: The integration contours in the 9 plane. The circles are poles at 6 = i7r/2 + 
sinh"i(27rnr/M). 

The correlation function of the disorder parameters can be presented in a similar fashion, 
(Mr,x)/i(0,0))= (10) 

°° rr2N+l 2Ar+l u.| i ( ) . 

cM^"^-''""^> E E n ^ n fefi . 

7v=o ^^^^ + gu-g2N+i i=i ^>j + / 

where now we have a sum over an odd number of g's. 

Eqs. (|7j) and (jlUj) are well placed to rewrite the correlation functions as an expansion in 
powers of exp[— /5M]. Focusing on the order-parameter correlation function, we rewrite Eq. 
((Tj) representing the A^-th terms in Eq. ((7j) as a sequence of contour integrals: 

1 r pT cosh ^i+la; sinh r pT cosh fijv+l^; sinh 6jv 

M y/^^ e/^cosh.._i -//^^ e/^cosh.._i ]lt^nl.^m-9,)/2], (11) 

with the contour C shown on Fig. ^ One can deform this contour into the superposition 
C Clow + Cup, where Ciow lies just above real axis while C^p corresponds to '^mO = tt — 0. 
As a result one obtains a sum of particle (lower contours) and anti-particle (upper contours) 
contributions: 

TV ^ „ n N-n 

n=l i=l j=l 

ntanh2[(^, _ 0^y2]l[tanh^m " 0'p)/2] x ncoth2[(0, _ _^ -^y^^^ ^^2) 

i>fc J>p i,p 

where /(e) = [e^^ — 1]~^ is the Bose distribution function and p = M sinh 6,e = M cosh 6. In 
Eq.lfT^. n and N — n are numbers of particles and antiparticles. Note that the singularities 



of the form factors are now removed from the real axis to one side of the contour, as was 



suggested in Ref. 



17[. A formula similar to Eq. (fT^ was proposed in Ref. |18j. However, 
there is a noticeable difference: the distribution functions turn out to be of the Bose type, 
although the quasi-particles of the Ising model are fermions. 

To obtain the casual Green's functions, one has to replace r with it. From now on we will 
be interested in the behavior of the Green's functions at longest real times t ^ M^^ and 
\x\ ^ M~^. Only those terms in Eq. (fT^ which contain both particles and antiparticles 
are singular. The leading contribution to the Green's function in each order comes from the 
vicinity of each singularity. Each 9in carries the thermal exponent e~^^ <^ 1. It is thus 
sufficient to restrict the summation in Eq. (jl2j) by terms with equal numbers of in and out 
rapidities, each in rapidity 6i being close to one of the out rapidities 6'j. There are n\ ways 
to choose a pair for each Of, this factor cancels one of the n\ terms from the denominator of 
Eq. ()12|). The integral over the difference 6i — 9'- can be evaluated with the resulting series 
summing up to the exponent, 

t)cr(0, 0)) = CM^/^e{t)e-^\''\<^^ exp[i?(t, x)], (13) 

where 

1 r o piMt(92-e|)/2+iMa;6»i2 



R(t,x>Q) = - I dd^de^e-^^'^^^'"'^^- 

TT^ J 



(012 + ioy 



TV 



This result reproduces the semi-classical computation found in Eq. (jH}. 

We now turn to how this result determines how we should proceed in the case of general 
integrable models. 



III. THE GENERAL CASE: TWO UNIVERSALITY CLASSES 

The form of Eq. (fT^ has a form intimately related to the double summed Lehmann 
expansion of Eq. 1. Apart the bosonic distribution function /(e) (unimportant in the limit 
T ^ M), Eq. (fT^ looks for all intent purpose as those subset of terms in Eq. 1 involving 
form factors (as given by Eq. 4) of exactly N particles or excitations. This then implies that 
the double summed Lehmann expansion is tractable. Or put another way, we can generalize 
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to correlation functions involving an operator O in other integrable models by using Eq. 
fjl2|) with the Ising model form factors replaced by the form factors of O. At temperatures 
much smaller than the gap, we replace the distribution functions by Boltzmann exponents. 
The asymptotics of the correlation function is then determined solely by the behavior of the 
form factors in the vicinity of the kinematic poles. As we have already pointed out, the 
existence of kinematic poles is a general consequence of the theory of form factors Q rather 
than a specific feature of the Ising model. These poles are due to annihilation processes 
and occur when rapidities of a particle and an anti-particle coincide. The behavior of form 
factors in the vicinity of the poles divide all integrable models into two universality classes. 
Accordingly the resolution of the problem of singularities should not be model-specific. 
Consider a generic matrix element of operator O: 



a\9i b) = {9i, . . . , 9n\0\6i, . . . , Om)ai,...a„;bi,...b^, (14) 



where 6i and 6j are rapidities and and bj are isotopic indices of in and out particles. One 
can express the form factor (|TH) in terms of canonical matrix elements. Let F^^. be the 
irreducible part of given by the connected diagrams. Only this par t is necessary for our 
calculations done to lowest order in exp (-M/3). According to Refs. ll3,Q,Q) K?r is equal 
to the matrix element between the vacuum and a certain excited state of n -|- m particles: 

F^{e,, a\e,, b),rr = {0\O\ei, . . . , -Itt + ^„ ± iO, ■ ■ ■ , -ivr + ± iO)^,,...,,^;,^,...^,; (15) 

where a indices are obtained from a ones by charge conjugation. The choice of whether the 
small imaginary part, iO, is positive or negative determines the exact form of the disconnected 
parts (here ignored) of the form factor in Eq. (fT^ . 

The axiom for the annihilation poles in the standard formulation is expressed in terms 
of the form factors of the type 

K,...,aAOu...,en) = {0\O\9n,...9l)a„,...,a,. (16) 

Namely, as was shown in Ref. 1^ (see also Ref. ^ in review), these form factors have poles 
at 9n — Oj = in, j < n, and the residues satisfy the following relation (see Fig. E}: 

iReSg^=0„_i+i,r-^f?,...,a„(^i'^2, . . . ,6'„) = :^5,...,a;_^(6'i, . . -,0^-2) (17) 



n— 3 

I _ ^2A^.,.^as<-^i'<(e^_^ - ^1) ( n s^'^iOn-i - 0^))s:zt'^:zl{en^l - ft„-2j 

i=2 
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FIG. 2: A graphical representation of Eq. H17|) written in terms of the irreducible form factors in 
Eq. (|14() (the rectangles). The up arrows represent rapidities 9i, while the down ones correspond 
to the rapidities — ivr + 6j. Circles at the intersections of the solid lines represent two-particle S 
matrices. 

Here S is the two-particle scattering matrix, C is the charge conjugation matrix, and <fa,o 
is the semi-locality index between the particle creation operator and O. We will define 
this index later. 

To demonstrate how Eq. ()17|1 can be used to extract the most singular behavior of the 
form factors, we combine it with (fT^ to obtain the following behavior of two- and three- 
particle form factors at the poles: 

ifl — e^'^^'^^-o) 

{e,a\0\ei,b) = — ^- ——^(0|O|0)5a,6 + disconnected part; (18) 

Oi — o 31 lO 

{ei,ai;e2,a2\0\e,a} = -i (0|C'|6'i - ivr, a^) 

17 — (72 ± lU 

A , (0^ - 0-,) - e*2'r<^>a,o A A , 
-1 j^^^ {OlOie, - ivr, a',) 

+ disconnected part. (19) 

Now it is appropriate to define the index (p. We will call two operators, Oi and O2, mutually 
non-local if the Euclidean correlator (. . . Oi{r)O2{0) . . .) is a multi- valued function of r. In 
particular, if under analytic continuation z — > ze'^'^^, z —>■ ze~'^'^^ [z = t + i x and z = t — i x), 
the correlator acquires only a phase, 27r(y9i 2, the two operators are called semi-local. The 
remarkable advantage of this definition is that the index (y^i 2 can be evaluated at small 
distances, where all correlation functions have a simple power law form. 

From Eq. ()18|1 it follows that an operator O with a finite vacuum expectation has 
annihilation poles only if ^a,o is not an integer, that is when O is non-local. For instance, a 
and /i in the Ising model are non-local with respect to the fermions (they acquire a phase, vr, 
after circling around the fermions in the r — x plane). However, if the operator can create an 

11 



odd number of particles, the annihilation poles may appear provided the S matrix is non- 
trivial. For instance, the vector field n in the 0(3) nonlinear sigma model belongs to the 

later category. Its form factors do have annihilation poles, although its vacuum expectation 
n 

value is zero 1]. 

We proceed with the calculation of the most singular parts of the form factors. We will do 
it explicitly for operators with (0|(9|0) 7^ 0. The maximal singularities come from the matrix 
elements which contain equal numbers of in and out particles. The contributions from all 
other matrix elements contain additional powers of either exp[— /JMl or exp[— itMl. Since 

f nfl 

form factors with different order of rapidities are related to one another (Refs. [11,12,121), 
it is sufficient to know the residues for one particular arrangement. For the form factors of 
the type in Eq. (fT^ . we focus on the residues at the poles at 9i = 6i. In this case one can 
use the crossing symmetry of the S-matrix and obtain from Eq. (fT7j) the relation for the 
residues of ()14j). This allows the convenient graphical representation of Fig. |21 

At this point, the computation divides itself into two cases. Each case marks out a 
separate universality class. We deal with each in turn. 



A. Ising-like Universality Class 

In the first case, Ising-like behavior results. This occurs if either of following conditions 
hold 

1. The S matrix is diagonal. 

2. The S matrix is not diagonal, but the semi-locality index, Lpa,o = 1/2, a = 1, . . . , a, is 
independent of the particle index, a. 

In these cases the residue is diagonal in the isotopic indices: = hi. For the diagonal S- 
matrix this is obvious. For the case if = 1/2 with a non-diagonal S-matrix, some additional 
reasoning is required. Namely, one needs to apply Eq. ()17p in successive steps. For the first 
step we obtain the residue in the two-particle form factor. This yields Eq. fll8|l diagonal in 
the isotopic indices. Then we use this result to calculate the coefficient at [(^1 — ^i)(02— ^2)]"^ 
in the four-particle form factor J^oi^^i — Itt ± iO, 62 — I'n ±: iO, 62, 6'i)a2ai62fei- follows from 
Eq. ()17|). the S-matrices enter into the expression for this coefficient in the vicinity of the 



12 



double pole in the combination 

Here we have used the crossing symmetry and unitarity of the S-matrices. As we see, this 
residue is again diagonal in the isotopic indices and equal to (—1)^(1 — e^'^^'^)^ = —4. 

Thus for these two cases, it is clear that a multiple application of formula (jl7j) yields the 
main singularity for general n as 

n 

. . . , . . . , e^) = {-ir - e'-'^--r U(0^ - 0P. ± i0)-^(0|O|0); 

p, i=i 

n 

= (-2i)"$^n(^.-^~P.±iO)^i(O|0|O), (21) 

p, i=i 

where the sum includes all permutations of {1, ...n}. The result for the correlation function 
of operators, O, with finite vacuum expectation values is then essentially the same as Eq. 
© but with 

a a 

exp[-(3e] ^sin^iTTipj) exp[-f3ej] = ^exp[-(3ej], (22) 
where the sum runs over isotopic indices of the excitations. 

B. Non-Ising Like Universality Class 

When both the S-matrix is non-diagonal and the semi-locality index, (j)a,o, is not simply 
±1/2, the behavior of the correlators is considerably more complex. Attendantly, the calcu- 
lations are more difficult due to the fact that (fa depends on the particle index. Nevertheless 
the maximal residue can be calculated, although only at small rapidities (in the given con- 
text, that is all we need). We take advantage of the fact that in integrable models with a 
non-diagonal S'-matrix, the S matrix has the universal limit S{6 — 0) = —P (the permuta- 
tion operator). This limit corresponds to a pure reflection. In general the crossover to pure 
reflectivity occurs at some model-dependent scale ^o- K may happen that this scale is much 
smaller than one; then universal dynamics occurs only for T <^ M6q < M when rapidities 
l^^l ^ ^0 determine the expressions for the correlation functions. At these temperatures, one 
can replace all S-matrices by permutation operators which leads to simplifications in Eq. 
dm) (see Fig. ^: 
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FIG. 3: A graphical representation of Eq. ()17|) for S=P. 

As we see from Fig. |Hl the isotopic indices of in and out particles must coincide. This 
allows us to simplify the notational representation of the form factor: 

Res,^^,-^^*'"^(^i, . . . , ai, ...aj^„, ...^i; a„, ...ai) = ...a„)(0|C>|0). 

Then the equation depicted in Fig. El becomes 



/(ai, ...an) = /(aa, ...a„) - e^'^'^"/(ai, ...a„_i) 



This recurrence relation can be solved explicitly with the result 



fc=l 



(l_ei<^^)(-l) exp (i ^ 

j=k+i 



n — I 
k-1 



where 0a = 2inpa,o- 

In this case the correlator, {O{x,t)O{0,0)), takes the form 



2tt 27r 27r 



((9(x,t)(9(0,0)) = {Orf:^j:\f{a,,...,an)H '^^^ ^^"^ ^^^^ 

n=0 ai 

/3M^cosh(^j- 



X exp 



X exp 



tMt ^(cosh(^j) - cosh(ej)) - iMx ^(sinh(^j) - sinh(^j)) 



n k 



■ on X] n //I 5 



n 



k=0 q= 



\ (e, -e,- w ^tk\i i(^P -^ + io) 



(23) 



(24) 



(25) 



where = (^). Again this expression arises either from using the form factors of O in 
conjunction with Eq. ()12|) or directly from the double summed Lehmann expansion (Eq. 

(I))- 



14 



Notice the particular arrangement of double poles in Eq. ()25p . As we have indicated 
the irreducible part of the form factor leaves certain ambiguity in whether to deform the 
poles into the upper or lower half plane. At present we have not been able to resolve this 
ambiguity. The advantage of this particular choice of poles is that the integral can be 
straightforwardly evaluated (under the same assumption as before that the most important 
contribution comes from small 9i and 9i). The result is 

(O(x,t)O(0,0)) = (0)'f;^^^El/(«i'- (26) 

71=0 ai 



S(x,t) = — e-*^^ / dOe-'^'^lx-tei. (27) 



where 

We now turn to the computation of the residues, /(ai, . . . , a„), of the poles. To simplify 
the discussion, we will consider only the case where the excitation spectrum of the theory 
constitutes a doublet: a = cr = ±1. This situation includes many important models, in 
particular the sine-Gordon model. In that case the semi-locality index of operator O with a 
soliton/anti-soliton equals ±0. We have, for instance, 

/(++) = (l-ei^)^ /(+-) = 2(1 -e-i*), /(-+)= 2(1 - e^*), /(— ) = (1 - e-i<^f . 

In this case Eq. reduces to 

J2 1/(^1, -^n)]' = 2"+^sin2(^) + 2sin\^)J2ctAClM-cos^r-^A . (28) 

cr,=±l I k>j ) 

To evaluate the sum, k < j, we rewrite it as follows: 



1 1 

^-^ COS(0) COSffl-^^ 



COS(0) 

For cos(0) > 0, we use the identities 



^ C'2n-2-r. (29) 



oo 



cos(0)l'=-^l = - / d^^^e*'^^'^-^^); 



V ^ = ^ ^^^'('^^ (30) 

^ 4(u; + m7r)2 + 52 8 sin2(cj) cos(0) + sin^(0/2) ' ^ ' 



m=—oo 
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where 5 = — In cos 0, to rewrite F as 

22n-4 _ Sin^(0) 

r = TTT / dcjsm (u;) — ^— -— , , , , , . (31) 

ncos{(f))jQ ^ sin2(cj)cos(0) + sm^((/)/2) ^ ^ 

If we finally use the identity to write the combinatorial factor, 6*2^2' ^ similar integral 
form 

= / d.;(cos(a;/2))2"-2, (32) 

we can, using Eq. ^I7\. represent the correlator as 

{O{x,t)O{0,0)) = - / dc.exp[-2g(x,t)sin^(u;)] . . (33) 

In Jq sm (cu) cos(0) + sm (0/2) 

It turns out that this relation holds even if cos(0) < 0. To establish this we must use instead 

the identity 

V - = - (34) 

„t^^ 4(cj + m7r)2 + 52 8cos2(cu)cos(0) + sin^(0/2)' ^ ^ 

where now 5 = — In | cos((/))|. We point out that as ^ tt, the expression for the asymptotics 
of the correlator (Eq. IHHj) reproduces the result for the Ising universality class. 

At a; = Eg. (1331) coincides with the correlation function obtained from the recent semi- 

nn 

classical analysis in Refs. jy, |21|. Both these papers explicitely assume that the two-particle 
S-matrix is equal to the permutation operator [2^. At finite x however, the semiclassical 
formulae are generically at variance with ours. In particular, at t = the semi-classical 
analysis rightly yields exponentially decaying correlation functions. This discrepancy pre- 
sumably originates from a different arrangement of poles in Eq. (j25|) . We hope to resolve 
this problem in a subsequent publication. As far as this paper is concerned, we will restrict 
the discussion to x = correlation functions. 



IV. LOW TEMPERATURE CHARGE DYNAMICS IN THE SINE-GORDON 
MODEL 

To provide an example of a physical problem which can be solved by the method developed 
in this paper, we consider the problem of charge counting statistics in the sine-Gordon model. 
The Lagrangian corresponding to this theory is 

^ = 7^(5,.$)' -"^'cos(7$). (35) 
16 



The cosine term is relevant at 7 < 1. The conserved sohton charge is defined as 

Q = ^J ^^9.,^- (36) 

Therefore Q{x,t) = (7/271) [$(2;, t) — $(a;,0)] is the charge that has passed through a point 
X during time t. One can calculate the probability distribution function of Q{x, t) using the 
results for the correlation function of operators exp[ir7$]. In this case (p± = ±(77/7). 

The sine-Gordon model is a good example to study as it encompasses both universality 
classes as 7 is varied. At the so-called refiectionless points 7^ = l/(n + l), n > 1, the 
sine-Gordon model model has a diagonal S matrix. At these points the model belongs to the 
same universality class as the Ising model. Away from these values of 7, the soliton-soliton 
scattering is not diagonal and the model falls into the second universality class (although for 
operators e^**"^/^ where the semi-locality index is ±1/2 we revert to the Ising universality 
class). In this case the physics is governed by the second universality class provided we are 
at a sufficiently low temperature. The rapidity scale governing this crossover is 

t'o^-cosh (cos( — ) + 2| sin(— )|), 

where p = 7^/ (1 —7^). If 7 is close to one of the refiectionless points, i.e. 7^ = l/{n + l) + 6, 
the crossover rapidity simplifies to 

n + 1 



ttS. 



n 



The corresponding temperature at which the crossover occurs is T* ~ ttM ^"^'^^-' 6. 

We now consider the counting statistics present in the sine-Gordon model when it falls 
into each of the two universality classes. 



A. Counting statistics at the refiectionless points 

At 7^ = l/{n + 1) we are in the Ising-like universality class. The correlation function 
then takes the form, as follows from Eq. (j^^ . 

^gi„${x,t)g-i^*(o,o)^ _ C2(7,r7)exp 

where e(p) is the soliton energy and the factor (7(7,77) is equal to the vacuum expectation 
value of exp[i?7$]. This formula remains valid even for 7^ 7^ (ra + 1)~^ provided r] = 7/2. 



-2 sin^ (7177/7) / — e~ 

J 



X ~ t 



de{p) 
dp 



(37) 
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The quantity C{'j,r]) was calculated in Ref. j23|: 

'aMT{l/{2 - 272))r((2 - 3^^)/{2 - 2^^))^^''' 



dy 

/o y 



4V? 
sinh^(2?77?/) 



X 



2 sinh(7^?/) sinh y cosh[(l — 7^)?/] 



2r/^e 



(38) 



Here we have assumed that the correlator, (^e^'^^^^'^^e satisfies the short distance 

normalization, 



^glr7*(a;,t)g-l»?$(0,0)^ 



4^2 



where a ^ is the effective cutoff in the theory. M ~ m^^^^ "^^^ is the soliton mass. In the 
field theory limit. Ma -C 1, and so the rj dependence in Eq. (IH8p is dominated by the term 
{Mafi . Thus the magnitude of the propagator has a strong maximum at small rj. 
We now examine the generating functional (exp[ici;Q(a;, t)]). This is given by 



) = (■ 



C(7,^)^exp[-2sin^(^) / %-^<^^\tdAp)\ 



(39) 



In the limit of large time the above integral finds its predominant contribution at small uj. 
We then can write 



(e^"'3(-,t)) = (aaM)(^)'exp 



cu'|t|/ro 



(40) 



where we have written C(7, |^)^ 
(T/7r)e-*'^/^. 

The Fourier transform of Eq. pn|) then yields a Gaussian distribution function for the 
charge, 



[aa 



M)(ir)' for some constant a and where I/tq 



P{Q) = {8{Q - Q)) = / dcue-i"'3(ei-Q(-.*)) 



~ a ^^'^ exp 



(41) 



where a = (t/ro) + ^ ln[l/(aaM)]. 
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B. Charge counting statistics at 7^ 7^ l/(n + 1) 



Here we consider the sine-Gordon model when it falls into the second universality class. 
As we have already mentioned, at 7^ 7^ + 1) there is a crossover temperature T* below 
which one can replace the S-matrix of slow solitons by the permutation operator. Then we 
may use Eq. finding that 

^gi.Q(x,i)^ ^ r du;exp[-AcuH/ro] , ^\ , . (42) 

The Fourier transform of this expression, after some manipulation, can be written as 

P{Q) = vrp/^ / dy , (43) 

J -00 ^/A+\y\ 

where i = t/ro and A = ^ ln(a~^/(a;M)). At small times, the integral sees its main 
contribution at small y, thus allowing us to write 

P(g),,.«^exp[-|(l-|j)], (44) 



In this case the distribution is that of a static Gaussian in Q with small time dependent 
corrections. At large times we can evaluate this integral in the saddle point approximation 
which yields 

/ 1 \ 1/6 r S / O \ 4/3-1 

piQm»i-{^) -p[-^(^) ]. (45) 

The reader may notice that quantum fluctuations play a significant role in modifying the 
Gaussian distribution of the counting statistics seen in the simpler Ising-like refiectionless 
case. 



V. CONCLUSIONS 



In this paper we have found universal expressions for the asymptotics of correlation func- 
tions of order parameter operators. These operators have T = nonzero vacuum expectation 
values thus breaking some discrete symmetry present in the theory. We restricted our con- 
sideration to theories with well defined asymptotic particle states (this excludes such models 
as the super sine-Gordon or the 0(2N-|-1) Gross-Neveu model). Among such models, we 
have found two characteristic types of order parameter dynamics. One type appears in the 
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models with either diagonal S-matrices (these include, among other models, the quantum 
Ising model and the sine-Gordon model at the refiectionless points) or for operators of non- 
diagonal theories with trivial semi-locality factors. These dynamics essentially coincide with 
the semi-classical dynamics described in ^|. Another type of dynamics is described by Eq. 
(|33|) and appears in the models where the zero momentum S-matrix corresponds to pure re- 
flection and the order parameter has a non-trivial semi-locality factor. Here the asymptotics 
of correlation functions, as determined by the form-factors, coincide with the semi-classical 



expressions found in Refs. [6| and [21| only at x = 0. We believe the origin of this difference 
is in the particular choice of the formfactor regularization. We hope to settle this problem 
in a subsequent publication. 
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